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A theoretical framework for analyzing stochastic data from single-particle tracking in complex or 
viscoelastic materials and under the influence of a trapping potential is presented. Starting from a 
generalized Langevin equation we found explicit expressions for the two-time dynamics of the tracer 
particle. The mean square displacement and the velocity autocorrelation function of the diffusing 
particle are given in terms of the time lag. In particular, we investigate the subdiffusive case. The 
exact solutions are discussed and the validity of usual approximations are examined. 
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1. INTRODUCTION 

The viscoelastic properties of complex fluids, like poly- 
mers, colloids and biological materials, can be derived 
from the dynamics of individual spherical particles em- 
bedded in it P, !3| . Particle tracking microrheology ex- 
periments [L, 2, 3;, ;4, 5, 6, 7, 8] are based on the obser- 
vation of the motion of individual tracer particles. In a 
typical microrheology experiment, particle positions are 
recorded in the form of a time sequence and information 
about the dynamics is essentially extracted by measur- 
ing the mean square displacement of the probe parti- 
cles d, im • Based on a generalized Langevin equa- 
tion with a memory function and assuming that inertial 
effects are negligible. Mason and coworkers [5, jj^i] have 
obtained a direct relation between the mean-square dis- 
placement of free tracer particles and the viscoelastic pa- 
rameters of the environment. It has been recently noted 
that the fluid inertia and the resulting memory effects 
become increasingly important when high-resolution ex- 
periments are performed [T^ . 

On the other hand, optical traps are increasingly used 
for position detection, with a wide range of applications 
in physics and biology 3, 13, 14, 15, 16, 17]. In an optical 
trap, the interaction between the laser and the trapped 
object can be approximated by a harmonic potential \18^ . 
However, the use of a trapping complicates the analysis 
of the obtained data, since the interactions with the vis- 
coelastic environment overlaps with the influence of the 
trapping force [J, [13]. For example, it has been noticed 
that neglecting memory effects leads to calibration errors 
of optical traps [T^ . 

It is now well established that when particles dif- 
fuse through a soft complex fluids or biological ma- 
terials, they exhibit anomalous subdiffusive behaviors 
[1, [m, [l^ [13, [llj. In this situation, the mean-square 



displacement exhibit a slow relaxation with the presence 
power-law decay in the range of large times. A theo- 
retical complete description of the behavior of a parti- 
cle in a complex medium and subjected to a harmonic 
potential can be formulated in terms of the generalized 
Langevin equation (GLE) containing a memory function 



13, S [13, [2^ • In a recent paper (23 , we have obtained 



analytical expressions for the evolution of mean values 
and variances in terms of Mittag-Leffler functions. How- 
ever, from an experimental point of view it is necessary 
to get expressions for the two-times correlation functions. 
For instance, the mean square displacement (MSD) can 
be expressed as 



p(r)= lim((X(i -I- r)-X(Of) 



(1) 



where \X{t + t) — X{t)\ is the particle displacement be- 
tween two time points, t denote the absolute time while 
T is the so-called lag time (26| . Alternative information 
about the experimentally observed diffusive behavior can 
be extracted from the normalized velocity autocorrela- 

H, 



tion function (VACF) [27|, defined as 



Cv{t) = lim 

t — '■OO 



{V{t + r)V{t)) 
{V^t)) 



(2) 
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Then, to calculate the MSD and the VACF one must 
know the behavior of the two-time correlations {X{t + 
T)X{t)) and {Vit + T)V{t)). In what follows we will in- 
vestigate the behavior of the MSD ^ and VACF ^ for 
harmonically bounded particle immersed in a viscoelas- 
tic environment. For this purpose, in Sec. 2 we present 
the corresponding generalized Langevin equation (GLE). 
The two-time dynamics is obtained, which enable us to 
calculate the MSD and VACF for arbitrary memory ker- 
nels. Section 3 is devoted to the study of the subdiffusive 
case. The analytical solutions are given and compared 
with the overdamped approximation. Finally, a Sum- 
mary of our results is presented in Sec. 4 . 
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2. DIFFUSION IN A HARMONIC WELL 
2.1. Formal solution for the GLE 

In what follows we consider the dynamics of a test 
particle of mass to, immersed in a complex or viscoelastic 
environment and simultaneously bounded in a harmonic 
potential well. The resulting motion can be described by 
the following GLE 

mXit)+m [ dt'-f{t-t')X{t')+mujlX ^ F{t), (3) 

"'0 

where ojq is the frequency of the trap, j(t) is the dissi- 
pative memory kernel, and the internal noise F(t) is a 
zero-centered and stationary random force with correla- 
tion function 



{F{t)F{t'))^C{\t^t'\)^C{r). 



(4) 



The integral term in ([3]) represents the dependence of 
the viscous force on the velocity history and the memory 
kernel ^{t) is related to the noise correlation function 
C{t) via the second fluctuation-dissipation theorem [2^ 



C{t) = kBTmj{t), 



(5) 



where T is the absolute temperature, and ks is the Boltz- 
mann constant. 

In what follows we consider the one-dimensional case, 
but our results can be easily extrapolated to the two or 
three dimensional case. The Langevin equation ([3]) can 
be formally solved by means of the Laplace transforma- 
tion. Taking into account the deterministic initial con- 
ditions xo = X{0) and vq = ^(0), the evolution of the 
Laplace transform of the position X{t) reads 



X{s) = Xq 



1 



ujI I{s) 



vo + -F{s)]G{s), (6) 

TO ' 



where F{s) is the Laplace transform of the noise. The 
relaxation function G{t) is the Laplace inversion of 



(7) 



where 7(5) is the Laplace transform of the damping ker- 
nel, and 



On the other hand, the Laplace transform of the ve- 
locity V{t) ~ X{t) satisfies that 



V{s) ^(vo + ^F{s) ] g{s) - xo uj'^ G{s) , (10) 



where 



g{s) = sG{s). 



(11) 



From Eqs. ^ and pU]) a formal expression for the 
displacement X{t) and the velocity V{t) can be written 
as 

Xit) = {Xit)) + - f dt'G{t-t')F{t'), (12) 



V{t) = {V{t)) + - / dt'g{t t')F{t') , (13) 



where 

(X(t)) = xo (1 - u^l I{t)) + «o G{t) , (14) 
{V{t)) = v^g{t)-xoulG{t), (15) 

are the position and velocity mean values evolution, re- 
spectively. 



2.2. Expressions for the MSB and VACF 

To calculate the two-time properties of the dynamical 
variables involved in the expressions of the MSD U) and 
VACF (121) we will make use of the double Laplace trans- 
form technique ^24*1. Then, from ^ and pH)) we have 

{X{s)X{s')) = xlx{s)x{s')+vlG{s)G{s') 
+ xoVo{x{s)G{s')+x{s')G{s)) 
+ \g{s)G{s'){F{s)F{s')), (16) 

{V{s)V{s')) = vlg{sW)+xW,G{s)G{s') 
- xoVQUj^{g{s)G{s') + g{s')G{s)) 

+ :^9{s)g{s'){F{s)F{s')), (17) 



(8) 



is the Laplace transform of 



where 



(18) 



I{t) = / dt'G{t') . 
Jo 



(9) is the Laplace transform of x{t) — ^ ^ ^o^it)- 



In the Appendix we show how the last term of Eqs. (fTB|) and (fT7|) can be calculated. Inserting expressions (|A.4p 
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and (jA.Sp into (fTBl) and ((T7)) and making a double Laplace inversion, we arrive at 

{X{t)X{t')) = xl x{t)x{t') + {vl - ^)G{t)G{t') + xo vo {x{t)G{t') + x{t')G{t)) 

+ ^ (7(0 + /(f) - I{\t - t'\)) - m)Iit') , (19) 

m m 

1 1 L 1 1 L 1 1 li 

- xovooj'o {9{t)G{t')+9{t')G{t)) (20) 
Finally, by considering time lags r > 0, from ([19]) and (pO|) we have 

((X(t + r)-X(t))') = ^M^/(r)-2xoi;oc^o'(G(i + r)-G(t))(/(t + T)-/(t)) 

m 

+ («o - ^) (G(t + r) - G(t))^ + (:,gc.2 - M) (/(t + ^) _ /(t))2 ^ (21) 

TO TO 

(n^ + r)T/(t)) = Mg(^) + („2_M:)g(t + ^)g(i)+^2(^2^2_M:)G(t + r)G(i) 

TO TO TO 



xo vo (g(i + T)G(t) + 5(t)G(t + r)) 



(22) 



Note that the analytical expressions ([2111 and ([22[1 are 
exact and valid for all absolute times t and time lags 
T. However, to evaluate the MSD ([J) and VACF ^ we 
must take the limit t —> oo. In this case, these expres- 
sions could be simplified as follows. Taking into account 
the usual assumption that the time-dependent frictional 
coefficient j{t) goes to zero when t — > cx) [IP] and using 
the final value theorem [30] one gets 



lim 7(t) = lim s'y(s) = . 

t — >oo s — >0 



(23) 



Noticing that the Laplace transform of the relaxation 
function I{t) defined through Eq. ([8]) is 



/(.) = 



(24) 



the application of the final value theorem and the use of 
condition ([23|) yields 



(25) 



I{oo) = I/wq, 
and using ^ and ([TT|) gives 

G(cx)) = g{oo) = 0. 



(26) 



Applying these conditions in order to take the limit 
i ^ cx) in ([2T|) and ([22| . and using the definitions ([T|) 
and ([2]) one finally obtain the simpler expressions 



TO 



and 



Gvir)=g{T). 



(27) 



(28) 



Taking into account ([^5]) and 
value of the MSD is given by 



p(oo) 



2kBT 



([^ . the equilibrium 



(29) 



while, as expected, the VACF decays to zero, i.e. 
Gv(oo) = 0. 

It is worth pointing out that in experimental realiza- 
tions the time lag is r^in < t < T,„ax, being T,„in the ac- 
quisition time interval and Tmax the measurement time. 
Moreover, if N is the number of steps n taken at inter- 
vals Tniin, only small values of rt (n < iV/10) are used. 
Therefore, it is important to obtain valid expressions for 
all observational time scales instead of getting only its 
behavior to large times. 

To conclude this section, we will find the extension 
for a trapped particle of the widely used Mason formula 
0, Q . Taking the Laplace transform of ([77]) and using 
the definition (|24p of the relaxation function I{t), one 
gets 



S7(s) 



2kBT 1 

TO Sp{s) 



(30) 



which gives a direct relation between the mean-square 
displacement of the particle and the memory kernel, from 
which the viscoelastic shear moduli of the medium can 
be obtained [l|. 



3. SUBDIFFUSIVE BEHAVIOR 

Notice that the previous results are valid for any mem- 
ory kernel that satisfy condition (|23p . On the other hand, 
it is well known that in the absence of active transport 
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the dynamics of the particle in a viscoelastic fluid or com- 
plex media is subdiffusive and thus the stochastic process 
presents a long-time tail noise. The most utilized model 
to reproduce a subdiffusive behavior is characterized by 
a noise correlation function exhibiting a power-law time 
decay [H, [sH El : 



C{t) = Cx 



r 



m-xy 



(31) 



where T{z) is the Gamma function [33]. The exponent 
A is taken as < A < 1 and the proportionality coeffi- 
cient Cx is independent of time but can depends on the 
exponent A. 

Using the ffuctuation-dissipation relation (O, the 
memory kernel j(t) can be written as 



lit) 



7A 



(32) 



r(l-A) ' 

where 7a = C\/kBT. Then, its Laplace transform reads 



l{s) = 7A s 



A-l 



(33) 



In this situation, the Laplace transform of the relax- 
ation function /(s) reads 



s + 7a s + 



(34) 



The complete temporal behavior of the relaxation func- 
tions I{t^G{t) and g{t) was previously obtained by us 
in Ref. [2^ . Using those results in ([7f)) and (|^ we have 



, , 2kBT ^ 

P(j) = > 

m ^ — ' 

fc=o 



(-1)^ 
fc! 



2k J2 



X E. 



(fc) 

2-A,3+Afc 



Cv{t) 



—rr- ('^or) 



fc=0 



fc! 



X 



(fc) 

2-A,l+Afe 



(-7Ar^-"), (36) 



where Ea^/s{y) is the generalized Mittag-LefSer function 
[3^ defined by the series expansion 



E, 



a, 13 



and E'^^'p (y) is the derivative of the Mittag-Leffler func- 
tion 



yj 



a>0, /3>0, (37) 



where the first term shows that the particle undergoes 
ballistic motion when time is very small .34] . The second 
term comes from the influence of the viscoelastic medium 
while the third term corresponds to the fact that the par- 
ticle begins to "see" the trap. The short times behavior 
of the VACF can be obtained in a similar way. In this 
case we get 



Cvir) 



7A 



-2-A 



r(3 - A) 



/ ,2 



(40) 



On the other hand, for "fxT^~^ > 1 the MSD and 
VACF can be obtained introducing the asymptotic be- 
havior of the Mittag-Leffler function [3^ , 



1 



y>0 



(41) 



yT{(3-a) 

into Eqs. (|55|) and After some calculations we have 



Pir) 
Cvir) 



2kBT 



1- Ex{ 



Jl (Pi 



7A 



^x^ 



(42) 

Exi-^r^), (43) 
t 7A 

denotes the one parameter 



where Ex{y) = Ex,v 
Mittag-Leffler function [33j . 

It is worth pointing out that these expressions can be 
also obtained discarding the inertial term in In 
this case we get 



/(.) = 



7A + ul 



„A-1 



'^o/7A 



(44) 



and using that the Laplace transform of the Mittag- 
Leffler function 33] 



Ea{-lt°')dt = 



+ 7 ' 



(45) 



one obtains expressions (|^^ and 

Finally, if » Ixl^'l the behavior of the MSD and 
VACF can be obtained using again the approximation 
(|4T|l . In this case we get 



Cvir) 



2kBT 



1-^ 



7A A(A 



1) 



r(i-A) 

.-(A+2) 



r(i-A) 



(46) 



(47) 



<%) 



dy 



ij + k)\y^ 



^.-j!r(a(j + fc) + /3) 



(38) 



Using the series expansions ([57]) and ([55)1 one can re- 
alize that the short times behavior of the MSD reads 



Pir) 



kuT 



27A 



.4-A 



r(5 - A) 



-^rH,(39) 
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showing a pure power law decay. 

In Figs. [1] and [2] we have plotted the MSD vs. time 
lag, using the exact solution (|35p and the approximations 
dH]) and (Uni)- Note that the exact solution exhibits a 
nonmonotonic approach to pioo), while the approxima- 
tions always present a monotonic behavior. Furthermore, 
even in the overdamped case the exact solution presents 
oscillations. The same differences can be observed in the 
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1.4 




2 4 6 8 10 



FIG. 1: (color online) MSD vs. time lag for A = 1/2, 7a = 1 
and ujo = 1.4. The solid line corresponds to the exact solution 
(|35|) , the dotted line to the approximate solution (|42|l and the 
dashed line to the asymptotic behavior (|46|l . 
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FIG. 2: (color online) Idem Fig. [T]for ljo = 0.8. 
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FIG. 3: (color online) Cv vs. time lag for A = 1/2, 7a = 1 
and uJo = 1.4. The solid line corresponds to the exact solution 
P6|) , the dotted line to the approximate solution (|43|l and the 
dashed line to the asymptotic behavior (|47|) . 
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FIG. 4: (color online) Idem Fig. [3] for ujq = 0.8 . 



behavior of Cy(r), as is evidenced in Figs. [3]and[31 In- 
terestingly, Burov and Barkai [ssj have recently arrived 
to similar conclusions examining the behavior of the po- 
sition correlation {X{t)X{0)). 

These behavior can be understood taking into account 
that the approximations (|42|) and (|43p only depend on the 
one parameter Mittag-Lcfflcr function E\{—uJqT''^/j\). 
On the other hand, it is known that the function E\{~t^) 
is a completely monotone function and tends to zero from 
above as t tends to infinity for < A < 1 [s^. Then, 
the approximate solutions are always monotonic for ev- 
ery value of luq and < A < 1. However, the exact 
solutions ([55)1 and are expressed as infinite sums of 

-^2-^ /3(~7-^ '''^ functions. In this case, the solutions 
can exhibit a nonmonotonic behavior as is displayed in 
the previous figures. 



4. SUMMARY 

In this work we have obtained the mean square dis- 
placement and the velocity autocorrelation function for a 
trapped particle and immersed in a complex or viscoelas- 
tic media. For this purpose, and starting from a suitable 
generalized Langevin equation, we have been able to de- 
rive analytic expressions for the two-times dynamics of 
the processes, valid for all absolute times and times lags. 
We have showed that the MSD and VACF can be ex- 
pressed as a simple expressions when the memory kernel 
goes to zero for large times. In particular, we have ex- 
amined the subdiffusive case, the which one is paradig- 
matic in the study of passive transport in viscoelastic 
media. In this case, exact expressions and valid for all 
time lags have been obtained in terms of Mittag-Leffler 
functions and its derivatives. The limit of short time lags 
are given in terms of the involved parameters. Finally, we 
have showed that the overdamped approximation, which 
means that the effects of inertia are neglected, can not 
reproduce the nonmonotonic dynamics present in the ex- 
act solutions. This result must be taken into account in 
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the analysis of the short and intermediate times dynam- 
ics where the MSD and VACF exhibit a relaxation plus 
an oscillatory behavior. 

In summary, we have presented a method to account 
for the effects of the trapping potential in the anoma- 
lous behavior of the mean square displacement and the 
normalized velocity autocorrelation function of a particle 
embedded in a complex or viscoelastic environment. 
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(^mit')) = Afi\t-t'\). 



(A.l) 



the corresponding double Laplace transform writes 



($(s)$(s')) =^ 



(A.2) 
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APPENDIX: 



Then, the Laplace domain version of the fluctuation- 
dissipation relation ([5]) reads ^2^] 



(F(s)F(s')) = kBTm 



7(s) +7(s') 



s + s' 



(A.3) 



To calculate the last term of Eqs. (|16p and P?|) we 
make use a relation given in Ref. [2J] . Given any station- 



After some algebra, and using the relations between 
the kernels I{t) G{t) and g{t) one can find that 



\ s' s s + s' 

g{s)g{s'){F{s)F{s')) = knTm ( ^^^1±4^ -g{s)g{s') - u;lG{s)G{s') 



ksTm (d{s)Gis')+u;lT{s)T{s')^ , (A.4) 

(A.5) 
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